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Abstract 

A micro-mechanics model for non-isotropic, open-celled foams is developed using an elongated 
tetrakaidecahedron (Kelvin model) as the repeating unit cell. The micro-mechanics model employs an 
elongated Kelvin model geometry which is more general than that employed by previous authors. 
Assuming the cell edges possess axial and bending rigidity, the mechanics of deformation of the 
elongated tetrakaidecahedron lead to a set of equations for the Young’s modulus, Poisson’s ratio and 
strength of the foam in the principal material directions. These equations are written as a function of the 
cell edge lengths and cross-section properties, the inclination angle and the strength and stiffness of the 
solid material. The model is applied to predict the strength and stiffness of several polymeric foams. Good 
agreement is observed between the model results and the experimental measurements. 

1. Introduction 

Previous studies on open and closed-cell foams have sought to establish a direct tie between the foam 
microstructure and the macro-level foam properties. Through careful consideration of the foam micro- 
structure and selection of a suitable representative repeating unit, equations for the foam density, elastic 
constants and strength have been written in terms of the micro-structural dimensions and the physical and 
mechanical properties of the solid material (Gent and Thomas (1959), Dement’ ev and Tarakanov (1970), 
and Huber and Gibson (1988)). 

To represent the foam micro-structure, many of these previous researchers used a tetrakaidecahedron, 
a fourteen-sided polyhedron comprised of six quadrilateral and eight hexagonal faces. The 
tetrakaidecahedron is widely known as the Kelvin foam model, as it was Thomson (1887) who, in his 
assessment of Plateau’s experiment, identified the tetrakaidecahedron (with slightly curved faces) as the 
only polyhedron that packs to fill space and minimize the surface area per unit volume (Gibson and 
Ashby (1997)). Zhu, et al. (1997), for example, adopted an equi-axed tetrakaidecahedron to develop 
equations for the foam Young’s modulus, shear modulus and Poisson’s ratio for isotropic, open-celled 
foam s . They assumed that the mechanical behavior of open-celled foams could be simulated by treating 
the edges of the cell faces as structural elements possessing axial, bending and torsional rigidity. 

Applying the principle of minimum potential energy to the deformation of the repeating unit, the 
equations for the foam elastic constants were written in terms of the cell edge length L, the edge cross- 
sectional area A, moment of inertia / and polar moment of inertia J and the Young’s modulus E and shear 
modulus G of the solid material. Using a similar set of assumptions, Warren and Kraynik (1997) 
developed similar equations for the Young’s modulus, bulk modulus and shear modulus for isotropic, 
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Figure 1. — Elongated tetrakaidecahedron repeating unit cell. 


open-celled foams. The more recent model by Gong, et al. (2005a) includes the effect of shear 
deformation and allows for the edge cross-sectional area to vary along the length of the edge. 

In many cases, the foam micro-structure is elongated in one of the three orthogonal directions causing 
the foam mechanical behavior to be non-isotropic. The micro-structure in closed-cell foams, for example, 
is often elongated in the rise direction due to the foaming and rising process. To treat non-isotropic 
foams, Dement’ev and Tarakanov (1970), Gong, et al. (2005a, b), Ridha, et al. (2006) and others have 
adopted an elongated tetrakaidecahedron (fig. 1) as the repeating unit cell, deriving equations for the 
elastic constants and strengths in the principal material directions. An elongated tetrakaidecahedron also 
packs to fill the space. It contains eight hexagonal faces, two horizontal square faces and four vertical 
diamond faces. The horizontal square faces have sides of length b and the diamond faces have sides of 
length L. The hexagonal faces have four sides with length L and two sides with length b. The inclination 
angle 0 defines the orientation of the hexagonal faces with respect to the rise direction as well as the 
obtuse angle of the vertical diamond faces. 

The size and shape of the elongated tetrakaidecahedron is uniquely defined by specifying the value of 
the three dimensions: b, L and 0. The above mentioned authors, however, have developed their equations 
for the elastic constants and compressive strengths of non-isotropic foams by imposing the restriction on 


the cell geometry that — = V2cos0 . This constraint forces the cell shape to be a function of the 


inclination angle only. Since, from a purely geometrical point of view, 0 and — may vary independently, 

we see no reason for this restriction on the cell geometry, other than to reduce the number of micro- 
structural measurements required to apply the equations and predict the foam behavior. As such, it is 
prudent to revisit the formulation of the previous authors and re-derive the equations for the elastic 
constants and strengths using the most general description of the elongated Kelvin model geometry. 

In this paper, we derive the equations for the Young’s modulus, Poisson’s ratio and strength for non- 
isotropic foams in the principal material directions. We follow closely the formulation from Zhu, et al. 
(1997), but adopt, as our repeating unit, an elongated Kelvin model with a geometry defined by three 
independent dimensions. Furthermore, we allow the edge cross-sections to assume any shape, but restrict 
our attention to edge cross-sections that do not vary along the edge length. In the application section, we 
make the simplifying approximation that the edge cross-sections are circular and apply the model to 
simulate the mechanical and strength behavior of several polymeric foams. In the final section, we extend 
the use of the model to closed-cell foams and predict the strength ratio of the five rigid polyurethane 
closed-cell foams studied by Huber and Gibson (1988). 
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2. Elongated Tetrakaidecahedron (Kelvin) Model 

2.1 Cell Aspect Ratio 

The size and shape of an elongated tetrakaidecahedron are uniquely defined by specifying the value 
of any three of the five dimensions L, b, 9, H, D, since the height H and width D of the unit cell is related 
to L, b, and 0 according to 


T/=4Tsin9 and D - 2Tcos9 + 42 b 


(1) 


The cell aspect ratio R = — is therefore 

D 


4Lsin0 

2Tcos9 + Jib 


( 2 ) 


There is a minimum value of 9, below which the unit cell in figure 1 is no longer elongated in the 

Z-direction. This minimum value of 9 is a function of the length ratio — , since as — becomes larger, the 

L L 

value of 9 must become larger in order for H > D and thus R > 1 . The equation for the minimum 9 in 
terms of the ratio — is derived in the Appendix. 


2.2 Foam Relative Density 

The relative density y is, by definition, the ratio of the foam density to the density of the solid 

P V 

material, y = — . The relative density may be written in terms of volumes as y = — , where V s is the 

P.s V 

volume occupied by solid matter and V is the total volume of the foam. Using the elongated 

2 

tetrakaidecahedron shown in figure 1 as a representative volume, the total volume is V—H D . The 
members that form the perimeter of the vertical diamond faces and those of length b that form the 
perimeter of the horizontal square faces are shared by the adjacent cells. Thus, they contribute only half 
their cross-sectional area to the repeating unit. All other members are completely contained within the 
boundaries of the unit cell. The volume of solid material is therefore V s = 1 6AL+8Ab , where, again, A is 
the edge cross-sectional area. Using the relations in equation (1), the relative density may be written 


y = 


2A(2L +b) 


Lsin9 


2Lcos9 + 


V2 bf 


( 3 ) 


2.3 Loading in the Y-direction (Perpendicular-to-Rise) 

We now seek to develop equations for the foam mechanical response and foam strength, for loading 
in the y (perpendicular-to-rise) direction, in terms of the micro-structural dimensions L, b, and 9, the edge 

cross-section area A and moment of inertia / and the modulus E and ultimate strength u lllt of the solid 
material. For this puipose, we establish the cartesian coordinate system shown in figure 2, where the 
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Z (rise direction) 


Z (rise direction) 



Figure 2. — Repeating unit cell for loading in the Y-direction (perpendicular to rise). 


Z-direction is oriented in the rise direction and the X and Y directions are in the plane perpendicular to 
the rise direction. We will develop the equations for the foam mechanical response and strength for 
loading in the Y -direction and note that due to symmetry of the repeating unit cell, the same set of 
equations will apply for loading in the X-direction. Furthermore, due to symmetry, we can use the unit 
cell shown in figure 2, which represents one-eighth of the tetrakaidecahedron shown in figure 1. 

We consider the deformation of this unit cell under the application of a uni-axial stress in the 
Y-direction ct vv which results in an extension in the Y-direction and the accompanying contractions in 

the X- and Z-directions. In order for the unit cell to be a representative repeating unit during deformation, 
we enforce the symmetry conditions on the member end point displacements: 

Up -Up- 0 Up — Uq — —u 

V C =v D = 0 v G = v H =v (4) 

w D -Wp-0 w B -w H --w 

where u, v, and w denote the displacements in the X-, Y-, and Z-directions, respectively, and u , v , and 
w represent the displacements of the unit cell at the unit cell boundaries. We also require that the 
deformation of the unit cell occurs with no rotation of the member end points. In addition, we note that 
due to similarity of members BC and FG and the similarity of members BF[ and DF, we have the 
additional conditions 


UjL[ — Up 
B _ 

V H -Vp 

w G = w B 


where vjj is the Y -direction relative displacement of point H with respect to point B, v H 


, B 

zV B + V H ■ 
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2 L sin 9 




* Y 


(a) Y-Z plane 


LcosO + b! V2 


(b) X-Y plane 


X 


Figure 3. — Unit cell deformation for loading in the Y-direction. 

The deformation of the unit cell is illustrated in figure 3. It is straightforward to write the strains in 
terms of the unit cell boundary displacements 


-2 u 


2v 


2Lcos9 + -Jib 
—w 


J yy 


2Z,cos9 +Jlb 


(5) 


2Tsin9 


Due to the loading on the unit cell and the restrictions on the displacements, the members CD and GH 
carry no axial load or bending moment and therefore they do not contribute to the strain energy. 
Furthermore, since the members BC and FG have the same length, orientation and stiffness and since the 
same is true for members BH and DF, the strain energy of the assemblage of members BC and BH are 
exactly equal to the strain energy of the assemblage of members DF and FG. Thus, the total strain energy 
of the deformation of the unit cell is 


U = l{U BC +U BH ) 


( 6 ) 


and we need only consider further the deformation of members BC and BH. 


Now since v H = v , we also have that 


A 


f B \ 

V B 

+ 

V H 

l V J 

V 

V J 


= 1 . Furthermore, it can be shown that, due to 


the relative stiffness of members BC and BH, the fractional displacements are 


v B _ 2Z.(l2/cos 2 9 + sin 2 9 ) 

v UlilLcos 2 Q + bj + a(2L 3 sin 2 0 + Z? 3 ) 

vj__ b(l2I + Ab 2 ) 

v \2l{lLcos 2 0 + &) + a(2L 3 sin 2 9 + Z? 3 ) 
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The strain energy for a single structural member carrying an axial force and bending moment is 
(Laible (1985)) 


U = 


N 2 L 


- + - 


2 AE 2EI J 0 


1 r L - 

— f M 2 (s)ds 
FT J 0 


( 8 ) 


where N is the axial force and M(s) is the equation for the bending moment written as a function of the 
local length coordinate s, where s is oriented along the member length. We note that all members in the 
one-eighth unit cell have a cross-sectional area of All, since they are on the boundary of the unit cell and 
are shared by two adjacent unit cells. Likewise, the moment of inertia of members BC and BH is 7/2, 
since, given the deformation of the unit cell specified by equation (4) and illustrated in figure 3, the axis 
of bending is perpendicular to the boundary surface for both members. Thus, the axial force and bending 
moment in members BC and BH may be written in terms of the member end point displacements as 

EA 

n bc =— ( v b cos0+w g sine) 


AT _ EA ( b\ 

* bh ‘uE [ uh * Vh > 


(9) 


M BC ( s )- 


f 3EI 6 El ' 


\ Zr L i 


(-VgsinO + vvgCosG) 


m bh( s )~ 


3-JlEl 3V2T7 


2ZT 


ir u H + v h ) 



/ \ 


( B \ 

Setting u H = -u and w B --w and making the substitutions v B -\ 

Vb 1 

1 R 

v and v H - 

vh 


l V ) 


V ^ y 


equation (9), the strain energy in members BC and BH are written in terms of the unit cell boundary 
displacements using equation (8) as 


U 


BC 


EA 
4 L 


f v B y 

vcos0- vvsinO 

2 

3EI 
+ ^ 


f V B^ 

vsin0 - vvcosO 

A v J 


I? 


[ V J 



EA 


( B 3 

vh 


2 

3EI 


f B 3 

vh 



— u + 

V 

+ 

u + 

V 

8 b 


y v j 


2 b 3 


y v ) 



(10) 


The work applied to the unit cell by the application of the applied stress a vv acting over the 
displacement v is 


W= 


a (2Lsin 0) 

yy v ' 


( 


Lcos 8+ 


V 


V2, 


( 11 ) 
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The product of the terms in square brackets in equation (11) represents the total force applied to the unit 
cell in the Y -direction, a force which is statically equivalent to the applied stress a vv acting over the area 


(2Lsin9) 


LcosO + —= 

S. 


Following Zhu, et al. (1997), we next apply the Theorem of Minimum Potential Energy (Love, 1944), 
which is 





dw 


Direct substitution of the expression for U from equations (6) and (10), however, leads to a trivial result. 
Since W - 2U and therefore U — W -U , we can write a mixed and more useful form of U which 
includes the applied stress and the unit cell boundary displacements. A more convenient form of the 
Minimum Potential Energy Theorem is then 

d{w-u) 0 d(w-u) 0 d(w-u) 0 (1; 

du dv dw 


Substituting equation (10) into equation (6) and substituting the result along with equation (11) into 
equation (12) leads to the following set of three simultaneous equations written in terms of the boundary 
displacements 


r \2EI EA A 

\ b 3 by 


u + 


12 El EA 

~b 


,3 


( V B ^ 

V H 


V V 2 


v = 0 


(13a) 


f . 




\2EI _ EA 

b 


V H 

V v J 


f . 


u + 


2 EA 2 A 24EI . 2 a 

cos'0-i- — — sm 0 


V v ^ 2 ^ 


+ 


V V J 


12 El EA 

h 3 b 


( B \ 2 

V H 


v v y 


2 EA 2AEI 


\ L 




V B 

cosOsinO 

A v J 



w = 2a vv Lsino(2Lcos9 + 42 b) 


(13b) 


7 EA 

12 El) 



/ 


— 


cosOsinO 

v + 

u 

L 3 J 

V l 7 2 


V 


12 El 2a EA ■ ^a 
— — cos 0 + — sm- 0 

V L 3 L 


w = 0 


(13c) 
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The axial force and maximum moment in the members of length L (members BC and FG) are 
obtained by substituting the expressions for v B and w B into the expressions for N BC and M BC in 
equation (9). The maximum moment occurs at the member ends, at either s = 0 or s = L. Recalling that 
-w B = w , the expression for v B and w B are obtained from the first expression in equation (7) and from 
equation (14). Substituting these into equation (9) yields 


M bc =M fg =^(Tsine) 2 (2Lcose + V26) 
N BC - N FG = g vv (2Lcos0 + V2Z?)zcos0sin0 


(18) 


The maximum stress in any member is cr max 


N M 

- — + — , where S is the foam edge section modulus. 
AS 


Substituting equation (18) leads to 


yy v , 2 ' 

LcosOsinG L sin" 9 

A 


2Lcos9 + 


42b 


(19) 


where S [ is the section modulus for the members of length L bending about the section neutral axis 
which is parallel to the unit cell X-direction. 

We assume that foam failure occurs when the foam stresses produce a maximum stress in any of the 
edges which is equal to the ultimate strength of the solid material, that is when tf max = o ll,t . We will 
assume that this criterion holds true whether the maximum stress is tensile or compressive. Therefore, the 
ultimate strength of the foam in the Y -direction, based on failure of the edges of length L, is given as 



_ ult 


L cos 9 sin 9 


- + - 


L 2 sin 2 9 

2 Sf 


2LcosQ + 42b 


The maximum moment and axial force in the members BF[ and DF is 


(29) 


M bh -M df 

N BH = N DF 


- — j=LbsinQ ( 2LcosQ + 42b ) 
2V 2 


= — ^Lsin9 (2Lcos9 + 42bj 
V 2 


( 21 ) 


This leads to 
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CT .vr =■ 


TsinO ZAsinO 
4l A 2V2 S b z 


2LcosO + 


42b] 


( 22 ) 


where S z. is the section modulus for members of length b bending about the section neutral axis which is 
parallel to the unit cell Z-direction. The ultimate foam strength, based on failure of these members of 
length b, is therefore 



_ult 


TsinO ZAsinO 
V2 A 2V25* 


2Zcos0 + Vzh] 


(23) 


In low density foams, bending stresses tend to be much more significant in contributing to material 
failure than the axial force contributions (Huber and Gibson (1988) and Ridha, et al. (2006)). Ignoring the 
axial terms in equations (20) and (23), the denominator in equation (20) will always be larger than the 

S ^ b 

denominator in equation (23) as long as the unit cell is elongated such that 4l sinO > — ) — . As a result, 

S h z L 

the edges with length L will fail at a lower applied stress a vv than the edges with length b and hence 
equation (20) will always yield a lower estimate of the ultimate strength than equation (23). 

2.4 Loading in the Z-direction (Rise) 

We now seek to develop an analogous set of equations for loading in the Z (rise) direction. For 
loading in the Z-direction, it is more convenient to use the repeating unit shown in figure 4. Since the 
members AC, BC, DC, and CE lie within the boundaries of the unit cell, and are not shared by an adjacent 
unit cell, it will be assumed that they have area A and moment of inertia Z. 



Figure 4. — Repeating unit cell for loading in the Z-direction (rise direction). 
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We define the unit cell displacements relative to point C, so that u c = v c = w c = 0 and impose the 
conditions on the displacements: 


U £) — ^ E — ^ 

v a =-v B =v (24) 

w A = w B = w 
w D =W E = -w 

Again, the deformation of the unit cell is assumed to occur with no rotation of the member end points. 

The displaced shape is illustrated in figure 5. Due to the anti-symmetry of the unit cell, u = v and 
therefore s xx = s vr . The relations between the unit cell strains and the boundary displacements are 


-2 v 

s xx — e vy — j 

2Z,cos9 + v26 


w 

TsinO 


(25) 


The members of length b are unstressed and therefore do not contribute to the total strain energy of 
deformation. Also, since the four members of length/, all have the same end displacements, they have the 
same axial force and bending moments. The total strain energy of deformation for the unit cell is therefore 
four times the strain energy of any of one these members, that is U — 4 U AC = 4 U BC — 4 U CD - 
The axial force and bending moment in member BC is 


EA 

N bc = —{v b cos9 + w g sin9) 


M bc ( s ) - 


'6EI \2EI ' 
v L 2 I? 


(~v B sin9 + w B cos9) 


(26) 



-JlL cos d+b 


Y 



i7/V 2 


2 L cos 6 + 42b 

(a) Y-Z plane (b) X-Y plane 

Figure 5. — Unit cell deformation for loading in the Z-direction. 
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Making the substitutions v B =-v and w B -w in equation (26) and using equation (8), the strain energy 
of deformation is obtained as 


U -4U bc = 2~^ (-vcos9 + w sin 9) 2 + (vsinQ + wcosQ) 2 (27) 

L L 3 

The work performed by the applied stress a zz acting over the displacement w on the top face and over 
the displacement -w on the bottom face is 


W - 2a zz (V2Lcos9 + &)~ w 
Applying the Minimum Potential Energy Theorem yields the set of equations 

IcosOsinG w - 0 


f T 12 1',') 

A cos - G + — — sin' G 

v - 

f 12/^| 

l L 2 ) 


l L 2 ) 


(28) 


(29a) 


f 2EA 24 El' 


L L 3 


cosGsinG v + 


24 El 2 2EA . 2 

— — cos 0 + sin G 

L 3 L 


w = <Uz(V 2T cosG + h) 


(29b) 


the solution to which yields 


v = -a. 


.{j2LcosQ + bj j — — 


2AE 24EI 


kosOsinO 


(30a) 


W - G, 


( nr , \2 f Esin 2 G 

: (V2LcosG + bj < 


L 3 cos 2 G 


2AE 24 El 


(30b) 


The Young’s modulus E z = 


f „ V 

< 3 ,, 


V S zz J 


is obtained by substituting the expression for w from 


equation (30b) into equation (25), which results in 


E z = 


24£7sinG 


2 „ 12/sin 2 G 

cos - G + 

[V2TcosG + b 

AL 2 



(31) 


£ yy 

Since , then u zx = u z>; . Using u zx = = — — , we get 


-8 


«zx = Uzy 




2vLsinG 


^2LcosG + s[2bj 


(32) 
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and upon substituting the expressions for the displacements from equation (30), equation (32) becomes 


V2 L 2 (aL 2 -12/)cos9sin 2 9 


u zx ^zy ]r 

12/Lsin Q + ALr cos" 9 V2Zcos0 + h 


Substituting the expressions from equation (30) into equation (26) yields for the axial force and 
maximum bending moment 

N bc = — ~ — — (V2Lcos9 + h) 2 


M BC 


a„ ZcosO 


(V2Lcos9 + h)' 


N M 

Since a max = — + — , we obtain 
A S 


2 A 4 Si 


sin9 Zcos9 /— r . , 

+ — V2Lcos9 + h 

O A a nl. 


and therefore the ultimate strength of the foam in the Z-direction is given as 


ult 

G zz — 


sin 9 Lcos9 

2^4 4 sf: 


V2Lcos9 + h 


2.5 Stiffness and Strength Ratios 


Defining the stiffness ratio as R E - — — , then equations (1), (2), (15), and (31) lead to 


2sin 2 9 + f — ) +^ 2 4|2 cos 2 9 + - 


L) AL 1 


i 12 1 1 

cos" 0 + — r-sin" 0 
AL 2 


Defining the strength ratio as R a - — and assuming that foam failure in the Y-direction occurs 

g I 11 ,., 


due to failure of the members of length L, then equations (1), (2), (19), and (35) lead to 
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( 38 ) 


2 Sy cos0 + JLsinO 

R o= R — L 

25 ^ sin0 + JLcosO 


2.6 Isotropic Foams 


In the special case of isotropic foams with an equi-axed repeating unit, we have the conditions H = D 
and E v = E z , or in terms of the aspect and stiffness ratios, R = Re = 1. Using equation (37), it is easily 

71 . 71 

shown that these two conditions are only satisfied when b = L and 0 = — . Setting b = L and 0 = — , both 

4 4 


equations (15) and (31) reduce to 


Ey=E z =- 


642EI 


r 4 

r 12 / 1 

L 

1 + y 


AL 2 _ 


= E 


(39) 


which is the expression obtained by Zhu, et al. (1997) for isotropic, open-celled foams. Also, the 
expressions for the Poisson’s ratios, equations (17) and (33), all reduce to 

* „ , [al 2 - 12 / 

u =0.5) — - — : — 

\AL 2 +12/ 

the expression for the Poisson’s ratio obtained by Zhu, et al. (1997) for the isotropic case. Furthermore, 

3 a 

the relative density becomes y = — — — — . 

2V2 L 2 


3. Application to Experimental Studies 

3.1 Relative Modulus Versus Relative Density in Isotropic Foams 

We first apply the equations derived in the previous section to simulate the variation of the relative 
E* 

modulus with the relative density y for isotropic foams. For isotropic foams, the relative modulus 

E 

follows from equation (39). Assuming that the edge cross-sections are circular with radius r, the relative 
modulus is 


* 

E 

E 


6y[l TZ 

4 


f r \ 4 




1 + 3 


( rY 


V-W 


(40) 
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and the relative density is 


Y = 


3rc 




\L) 

2V2 


(41) 


Y 

Solving equation (41) for the ratio — and substituting the result into equation (40) leads to 

L 

E* _ 4V2 y 2 
E 3 n + 2y[2y 


(42) 


By comparison, the relative modulus for a square cross-section is — r r , which is similar to 

3 (3 + 2V2 yj 

equation (42), except that n is replaced by 3. 

The variation of the relative modulus with the relative density given by equation (42) is plotted in 
figure 6 along with the experimental measurements reported by a number of previous researchers. The 
results are plotted on a log-log scale. Equation (42) provides a good correlation with the experimental 
results, particularly at low relative densities. 



Figure 6. — Relative modulus plotted versus relative density for isotropic, 
open-celled foam. Results from previous experimental studies are 
plotted along with the relation given by equation (42). 
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3.2 Stiffness and Strength Ratio Versus Cell Aspect Ratio and Relative Density for 
Non-isotropic Foams With a Restricted Unit Cell Geometry 

The stiffness ratio and its dependence on the cell aspect ratio has been measured and reported by 
previous researchers for several non-isotropic polymer foams and the variation of the strength ratio with 
cell aspect ratio for Porolon was measured by Polyakov and Tarakanov (1967). In this section, we seek to 

rewrite the equations derived in section 2 for the ratios Re and R a , and obtain expressions for these ratios 
in terms of the two variables R and y, so that we may compare the model predictions with the 
measurements made by previous researchers. 

In view of equations (37), (38), and (2), it is clear that the modulus ratio Re and the strength ratio R a 
are functions of 0, b, L and the cross-section properties. In order to rewrite the expressions for Re and R a 
in terms of R and y, it is necessary to impose an additional condition on the unit cell geometry to reduce 
the number of u nk nown micro-structural dimensions by one. We impose the condition on the elongated 

Kelvin model that — = V2 cosO , which was originally suggested by Dement’ev and Tarakanov (1970) 

L 

and also assumed by Gong, et al. (2005a). This constraint leads to D = 4TcosO and R = tan 0 and thus 


cos# = 


sin 9 = 


Furthermore, the relative density, equation (3), now becomes 

2(1 + 7? 2 ) 3/2 +yf2^ + R 2 ) 


Substituting the relations in equation (43) into equations (37) and (38) yields 


2 R 2 + 


-^— + ^■[2 + ^1 + R 2 

KLL 2 


1 + ™R 2 


R a =R 


+R 

AL 


( 2 

R + 1 

AL 

v y 



NASA/TM— 2007-214931 


16 




_127__ JVj 2 24Ky 

al1 ^ 27r(l + ^ 2 ) 3/2 +V2 Tt(l + 7? 2 ) 



Substituting the relations in equation (47) into equations (45) and (46) leads to 



We note, once again, that an expression similar to equation (48) is obtained if a square cross-section 
is assumed, except, in that case, n in the numerator and denominator is replaced by a 3. 

Equations (48) and (49) are analogous to equation (1) and (4) in Huber and Gibson (1988). Although 
considerably more cumbersome, equations (48) and (49) are more inclusive than Huber and Gibson’s 
relations, as they account for both axial and bending deformations of the cell edges and they include the 
effect of the relative density on the stiffness and strength ratios. 

Using equation (48), the modulus ratio is plotted versus cell aspect ratio and relative density in 
figure 7. The experimental results reported by a number of previous researchers are also included. Note 
that the stiffness ratio is a weak function of the relative density, particularly as the cell aspect ratio 
approaches unity. As the cell aspect ratio approaches unity, all the curves must collapse to the same point, 
that is, at R - 1 and Re = 1 . 

The variation of the strength ratio with the cell aspect ratio, given by equation (49), is plotted in 
figure 8. The variation of the strength ratio with the cell aspect ratio for Porolon as reported by 
Dement’ ev and Tarakanov (1970) (originally reported by Polyakov and Tarakanov (1967)) is also plotted 
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Cell Aspect Ratio, R 

Figure 7. — Modulus ratio plotted versus cell aspect ratio and relative 
density. Results from previous experimental studies are plotted 
along with that given by equation (48). 



Figure 8. — Strength ratio plotted versus cell aspect ratio. Results on Porolon 
from Polyakov and Tarakanov (1967) as reported by Dement’ev and 
Tarakanov (1970) plotted along with equation (49) fory = 3 percent. 

in figure 8. A relative density of 0.03 was used to plot equation (49), as this was the relative density 
reported in the experimental study. Clearly, equation (49) provides a reasonably-close match with the 
measured results, given the small amount of the experimental data available. 

3.3. Analysis of Polymer Foams Using a General Description of the Unit Cell Geometry 

We now apply the equations derived in section 2 to simulate the strength behavior of the rigid 
polyurethane foams reported by Huber and Gibson (1988). We will allow for the most general description 
of the unit cell geometry, which requires us to specify three dimensions to describe the unit cell geometry 
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plus a description of the edge cross-section. We will assume a circular edge cross-section and therefore 
we need only to define the cross-section radius. Using the values for the cell width, the cell aspect ratio, 
the relative density and the modulus ratio reported in Huber and Gibson (1988), we use the equations 
derived in section 2 to calculate the unit cell dimensions and the edge cross-section radius for each of the 
five rigid polyurethane foams. Using these four dimensions, the strength ratio for each of the five foams 
can be calculated using a modified form of equation (38). 

Now since the rigid polyurethane foams studied in Huber and Gibson (1988) were closed-cell foams, 
their microstructure includes cell faces which are stressed and deform with the application of the applied 
stress. As such, the faces also contribute to the foam mechanical behavior. In addition, the edge cross- 
sections are not really circular, but rather they resemble a three-cusp hypocycloid (Gong, et al. (2005a)). 
As a result, the edge radii that are calculated for the foams are really an equivalent cross-section radius. 
That is they represent the radius of a circular cross-section possessing a structural rigidity that is 
equivalent to the combined structural rigidity of the cell faces and the three-cusp hypocycloid edges. 

From equation (2), we have 


b 4sin0-2/?cos0 

L JlR 


(50) 


JJ 

Using D - — and equation (1), equation (3) can be written as 
R 


AR- 


2 + - 


8 L 2 sin 3 0 


(51) 


Assuming a circular edge cross-section with radius r and upon substituting equation (50), equation (5 1 ) 
can be rewritten and rearranged as 






4V2ysin 3 Q 

hr[42R + 2 sin 9-/? cos 0 


(52) 


12 / 

Also, assuming a circular edge cross-section allows us to replace — — with 3 


/ r 5 2 


AL Z 


and rewrite 




equation (37) as 


6R 2 cos 2 0-12/?£ sin 2 0 


M 

z 

1 ~ ^2 1 


M 

1 

\ T~\ 2 1 

( b ) 

Uv 

+ 3R Z \ 


,LJ 

+ R 



+ 2/? 2 sin 2 Q-4R e cos 2 0 = 0 


(53) 


In order to solve for L, b, 0, and r, we could substitute equations (50) and (52) into equation (53) and, 
using the identity relation cos' 0 = 1 - sin' 0 , obtain an explicit algebraic equation in sin 0 which can be 
solved to obtain the value of 0. The resulting equation is, however, quite cumbersome. Instead, we choose 
to solve the set of equations (50), (52), and (53) in an iterative manner, whereby these equations are 
solved in series within each iteration. 
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The steps in the solution approach were to make an initial guess at the value of the inclination angle 

b r 

0, substitute this value into (50) and (52) and solve (50) and (52) for — and — , respectively. The values 

L L 

b r 

of — and — were then substituted into equation (53) and equation (53) was solved for 0. This process 
L L 

was repeated until the difference in the value of 0 between two successive iterations was within an 
acceptable tolerance. The standard bisection method (Press, et al., 1992) was implemented to solve 
equation (53) numerically for 0. The method converged usually within thirty iterations given an initial 
bracketed guess for the inclination angle between 40 and 15°. Once the final value of 0 is obtained from 

RD 

the iterative solution, the value of L is determined from L — and the value of b and r is calculated 

4sin0 

from equations (50) and (52), respectively. 

The cell width, relative density, cell aspect ratio, modulus ratio and the measured compressive 
strength ratio for each of the five rigid polyurethane foams reported in Huber and Gibson (1988) are listed 
in table I. We have also listed the equivalent edge cross-section radius, the inclination angle and the 
lengths of the edges L and b which were obtained from each iterative solution. Notice that the equivalent 
cross-section radii seem to correlate with the measured relative densities, and the inclination angles seem 
to correlate with the measure modulus ratios. Also, although the elongated tetrakaidecahedron shown in 
figure 1 is drawn so as to imply that the edges of length b are shorter than the edges of length L, this is not 
strictly true. Indeed, one of the foams listed in table I has b = L and another has b> L. 


TABLE I.— RESULTS FROM THE EXPERIMENTAL STUDY ON RIGID POLYURETHANE FOAMS BY 
HUBER AND GIBSON (1988) ALONG WITH THE RESULTS FROM THE NUMERICAL SOLUTION 


Experimental data for rigid polyurethane foams 
from Huber and Gibson (1988) 

Solution results 

Cell width, 
D, 
mm 

Relative 

density, 

Y 

Cell aspect 
ratio, 

R 

Modulus 

ratio. 

Re 

Exp. comp, 
strength 
ratio, 3 

Ra 

Edge cross- 
sectional 
radius, 

r, 

mm 

Inclination 

angle, 

0, 

deg 

Edge 

length, 

L, 

mm 

Edge 

length, 

b, 

mm 

Predicted 
strength ratio, 

R a 

0.20 

0.027 

1.405 

2.518 

1.730 

0.0070 

53.91 

0.087 

0.069 

1.880 

0.20 

0.053 

1.470 

3.312 

1.943 

0.0100 

56.63 

0.088 

0.073 

2.132 

0.20 

0.080 

1.230 

2.256 

1.484 

0.0117 

53.29 

0.077 

0.077 

1.582 

0.20 

0.107 

1.275 

1.475 

1.152 

0.0137 

49.26 

0.084 

0.064 

1.447 

0.20 

0.133 

1.190 

2.122 

1.455 

0.0149 

52.98 

0.075 

0.078 

1.498 


“Experimental compressive strength ratios were calculated from the plastic collapse stresses reported in Huber and 
Gibson (1988). 


Assuming a circular cross-section of radius r, equation (38) becomes 

| + tan0 


f r \ 


\2Lj 


R ° R r „ v 


2 L 


tan0 + l 




(54) 
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2.7% 5.3% 8.0% 10.7% 13.3% 


Relative Density, y 

Figure 9. — Comparison of the measured and predicted strength ratio 
plotted versus the foam relative density for rigid polyurethane foams. 
Experimental data from Huber and Gibson (1988). 


The strength ratio was calculated for each foam using equation (54) and the values of r, 9, and L that were 
obtained from the iterative solutions along with the measured values of R. The predicted strength ratios 
are listed in table I; a comparison of the measured and predicted strength ratios is shown in figure 9. The 
predicted strength ratios are within 10 percent of the measured strength ratios for all cases except for the 
10.7 percent relative density foam where the error is close to 26 percent. It should be mentioned that the 
value of the cell aspect ratio reported by Huber and Gibson (1988), for this foam, was brought into 
question, as they reported difficulties with the microscope while measuring the micro-structure of this 
foam. Further, the cell aspect ratio for this foam does not appear to be consistent with the measured 
strength and stiffness ratios. The measured strength and stiffness ratios for the 10.7 percent relative 
density foam were the lowest of the five foams, whereas the cell aspect ratio was the median of the five 
values. Ignoring the results for the 10.7 percent relative density foam, we can conclude that the equations 
derived in section 2 as well as the iterative solution were successful in predicting the compressive strength 
ratios for the closed-cell foams studied by Huber and Gibson (1988). 

Finally, we note that, for all five cases, the value of the ratio is not equal to the value of V2 cosO , 

reinforcing our notion that the restriction on the unit cell geometry used by the previous researchers is 
unfounded and adopted merely for the sake of convenience. 


4. Concluding Remarks 

The formulation by Zhu, et al. (1997) has been revised to include the mechanical and strength 
behavior of non-isotropic foams. Equations for the foam Young’s modulus, Poisson’s ratio and strength 
in the principal material directions were obtained by adopting an elongated Kelvin model as the repeating 
unit cell. These equations were written in terms of the edge lengths and edge cross-section properties, 
the inclination angle and the strength and stiffness of the solid material. The micro-mechanics model 
is developed from the most general description of an elongated Kelvin model, as it requires three 
independent dimensions to describe the unit cell geometry as well as a description of the edge 
cross-section. 
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The model was applied to simulate the variation of the relative modulus with relative density in 
isotropic foams as well as the variation of the modulus ratio and strength ratio with cell aspect ratio in 
non-isotropic foams. In all cases, the model results were in good agreement with the experimental 
measurements. The model was also applied to simulate the strength ratio in closed-cell rigid polyurethane 
foams. Here, also, the model results were in good agreement with the measurements, as the predicted 
strength ratio was within 10 percent of the measured strength ratios for all but one of the five foams. 

In closing, it is worth noting that by adopting an elongated Kelvin model with the most general 
geometry, a more detailed description of the foam microstructure is required in order to apply the 
resulting equations and predict the foam behavior. More specifically, it is now necessary to obtain four 
separate physical and mechanical measurements of the foam in order to apply the equations. Aside from 
this added burden, the model is an improvement over the previous models, since it is capable of more 
closely representing the foam micro-structure for a wider range of foam materials. Thus, the resulting 
equations should more accurately simulate the foam behavior for a wider range of foams. 
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Appendix 

Relation Between the Minimum Allowable Inclination Angle and 
the Length Ratio bIL for an Elongated Tetrakaidecahedron 

If the height H and width D of the tetrakaidecahedron are equal, the expression in equation (1) leads 
to 


2sin0-cos9 = 


4il 


(Al) 


For an elongated tetrakaidecahedron ( H > D), we have the inequality 

b 


2sin0-cosO > 


<JlL 


(A2) 


Using the trigonometric identity, cos 9 = Vl-sin 2 0 , equation (Al) can be rearranged and rewritten as a 
second-order polynomial in sin 0, that is, 


5sin 2 0-2 a/ 2— sin0 + 
L 




y 2L- , 


= 0 


(A3) 


The solution to equation (A3) is obtained using the quadratic formula resulting in the two solutions 

(A4a) 


. . 4ib V 2 I b 2 

sin0i = + , 10 — 


sin9? = 


5 L 10 


V2 b V2 


10 — —rr 


5 L 10 V L 2 


(A4b) 


where both roots are real provided — < \/To and where the two are identical when — = Vl() . The 

L L 


solutions listed in equation (A4) are plotted in figure Al. 

It is easily shown that equation (A4a) is the solution to equation (Al) and that equation (A4b) is the 
solution to 


2sin0 = -cos0+ L- (A5) 

V2 L 

Since equation ( A5) has no physical significance here, we will ignore the latter of the two solutions in 
equation (A4). 
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The plot of equation (A4a) in figure A1 defines the value of sin 0 as a function of the length ratio — 

for any tetrakaidecahedron with H = D. As such, it defines the lower bound on the inclination angle for all 
possible elongated tetrakaidecahedron, since sin 0 for any elongated tetrakaidecahedron must lie above 
the upper curve in figure Al. 

2 b 

We note that, in the range — -= < sin0 < 1 , equation (A4a) yields two possible values for the ratio — 

V5 L 

for each value of sin 0. The values of — > 2 -Jl , however, violate equation (A2) since 


2>2sin0-cos0 for all 0< 


71 


Hence, for an elongated tetrakaidecahedron with 0 < — , the length ratio — must be less than 2^2 . 

2 L 

The lower bound on the inclination angle is therefore given by 


0 > Arc sin 




+—J10 — - 


5 L 10 


(A6) 


which is valid over the domain 0 < — < 2 V 2 . 

L 
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